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Abstract 

We consider a new partial differential equation, of a similar form to 
the Camassa-Holm shallow water wave equation, which was recently 
obtained by Degasperis and Procesi using the method of asymptotic 
integrability. We prove the exact integrability of the new equation by 
constructing its Lax pair, and we explain its connection with a negative 
flow in the Kaup-Kupershmidt hierarchy via a reciprocal transforma- 
tion. The infinite sequence of conserved quantities is derived together 
with a proposed bi-Hamiltonian structure. The equation admits exact 
solutions in the form of a superposition of multi-peakons, and we de- 
scribe the integrable finite-dimensional peakon dynamics and compare 
it with the analogous results for Camassa-Holm peakons. 



1 Introduction 

Since its discovery [fl| there has been a considerable amount of interest in 
the Camassa-Holm shallow-water equation, 

u t + 2ku x + ju xxx - u xxt + 3uu x = 2u x u xx + uu xxx . (1.1) 

It was originally derived as an approximation to the incompressible Euler 
equations, and found to be completely integrable with a Lax pair and as- 
sociated bi-Hamiltonian structure 0. It was subsequently shown to have 
a hodograph link to the KdV hierarchy and an interpretation within the 
framework of hereditary recursion operators ||. In the CH equation with 
linear dispersion was shown to be one full order more accurate in asymptotic 
approximation beyond Korteweg-de Vries (KdV) for shallow water waves. 
Yet, it still preserves KdV's soliton properties such as complete integrability 
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via the inverse scattering transform (1ST) method. 2+1-dimensional ana- 
logues have also been studied |L4|, jig] . A particularly unusual feature of 
(1.1) is that in the dispersionless limit k — > it admits peaked solitons or 



peakons. The single peakon takes the form 

u(x,t) = ce-\ x ~ ct \, (1.2) 

while the iV-peakon solution is just a simple superposition, 

JV 

u(x,t) = J2Pj(t)e- lx - qM , (1-3) 

3=1 

where the canonical positions qj and momenta pj satisfy a completely inte- 
grable finite-dimensional Hamiltonian system. 

The current work was motivated by the question of which equations of 
similar form to ( |1.1| ) are integrable. In the recent study the method of 
asymptotic integrability was applied to the family of third order dispersive 
PDE conservation laws, 

ut + c u x + ^u xxx - a 2 u xxt = (ciu 2 + c 2 u 2 x + c 3 uu xx ) x , (1.4) 

whose right-hand side is the derivative of a quadratic differential polynomial. 
Within this family, only three equations that satisfy asymptotic integrability 
conditions up to third order are singled out, namely KdV (a = c 2 = C3 = 0), 
Camassa-Holm (ci = — ^cs/a 2 ,c 2 = C3/2) and one new equation (c\ = 
— 2cs/a 2 , c 2 = C3) which can be scaled to the following form: 

2 ( 9 3 \ , . 

u t + u x + 6uu x + u xxx - a I u xxt + -u x u xx + -uu xxx 1=0. (1.5) 

(Unfortunately, the equation in [B| appears with the typographical error 
3/2 — > 2/3.) However, by its nature the method of Jl only isolates nec- 
essary conditions for integrability and is not sufficient to prove it. In the 
following we demonstrate that the new equation ( |1.5| ) is integrable by ex- 
plicitly constructing its Lax pair. 

A simple observation that proved very useful in our investigations was 
that by rescaling, shifting the dependent variable and applying a Galilean 



boost, the equation (1.5) may be transformed to the dispersionless form 



ut - u xxt + 4:uu x = 3u x ii xx + uu xxx . (1.6) 

but integrability is unaffected by such transformations. (We have set the 
parameter a = 1, but it will sometimes be convenient to reintroduce it when 
considering various limits.) The next key observation was that in this form 
the equation (|0]) also has the single peakon ( |1.2| ) as an exact travelling 
wave solution. This led us to consider the family of equations 

u t - u xxt + {b + l)uu x = bu x u xx + uu xxx . (1.7) 
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for real parameter 6, which includes both Camassa-Holm (6 = 2) and the new 



equation (L6) (6 = 3) as special cases. It turns out that all the equations 



in this family possess not just the peakon solution (|1.2|) but also the multi 



peakon solutions (1.3). In the case of arbitrary b the pj, qj are not canonical 



variables, but satisfy the dynamical system 



tu r\ . 9G N 



where the generating function Gjy is 

1 N 

j,k=l 

this takes the canonical Hamiltonian form only in the special case b = 2 
(Camassa-Holm) . 

The asymptotic analysis in [B implies that 6 = 2,3 should be the only 
possible integrable cases within the family ( |1.7| ). In a class of PDEs 
including all of the equations ( |1.7|) was tested for integrability using Painleve 
analysis. However, because both Camassa-Holm equation ( |1.1| ) and the new 
equation (|0]) are examples of integrable systems with algebraic branching 
in their solutions (the weak Painleve property of ]18|), they were explicitly 
excluded by the various Painleve tests applied in [1(|, and in fact all of 
the equations in that class failed the combination of tests. The authors of 
JlO[ noted that the (strong) Painleve property is destroyed by changes of 
variables, and thus a transformation may be required before applying the 
test. For the case of Camassa-Holm a reciprocal transformation provides a 
link with a negative flow of KdV (see || [l2|]), which restores the Painleve 
property. In a forthcoming article S we shall generalize this transformation 



for the whole family (1/7). Applying Painleve analysis to the transformed 
family of equations then isolates the cases b = 2,3 as the only candidates 
for integrability. 

In this note we restrict our attention to the new integrable equation 



(1.6). We first present the associated linear system with a third order Lax 
operator. We explain the connection with a negative flow of the Kaup- 
Kupershmidt hierarchy pi] , and derive infinitely many conservation laws 
and a proposed bi-Hamiltonian structure. The dynamics of multi-peakon 
solutions, described by the finite-dimensional system ( |1.S| ) in the particular 
case 6 = 3, is also analysed. 

2 Reciprocal transformation 

We consider the equation ( |1.6[ ) which we write in the form 

mt + rn x u + 3mu x = 0, m = Cu := (1 — d%)u. (2-1) 
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In this form it is clear that we have the conservation law 

{m l % = -{m l '\) x . (2.2) 

The construction of the other conservation laws (infinitely many) will be 
presented in Section 4. 

Our strategy is to make use of a reciprocal transformation which maps 
between systems of conservation laws [15, [L9], pOfl , since it is known that such 
transformations can be used to connect systems with algebraic branching 
in their solutions to equations which have only pole type singularities. The 
relevant example here is the Camassa-Holm equation, which has a reciprocal 
transformation to the first negative flow in the KdV hierarchy, namely the 
equation 

R ■ V T = 0, (2.3) 

where R = d\ + AV + 2Vxd x l is the KdV recursion operator (for more de- 
tails see [12] and references). Camassa-Holm has weak Painleve expansions 
with cube root branching [p3], while KdV solutions have only double pole 
singularities. 

To define a reciprocal transformation for the equation (2.1) we introduce 
the dependent variable 

p = - m 1/3 (2.4) 

which is related to u by 

p 3 = (d 2 x -l)u. (2.5) 



We then use the conservation law ( [2. 2D in the form 

p t = ~(pu) x (2.6) 
to define the new independent variables X, T via 

dX = pdx — pudt, dT = dt. (2-7) 
By transforming the derivatives we find the new conservation law 

(p-!) T = u x . (2.8) 



Replacing d x by pdx in ( |2.5[ ) and using ( |2.8| ) to eliminate derivatives of u 
leads to the identity u = —p(}ogp)xT —p 3 - Hence, (|2.8|) can be written in 
terms of p alone, viz: 

{p^)t + (p(logp)xT+p 3 )x = 0. (2.9) 

It then turns out that this equation for p can be written in conservation 
form in another way, namely 



Pxx . P X 1 \ 3 ^ 2 A /r> 1n s 



p 2p 2 2p 2 J T 2 
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The equation ( |2.10| ) is most conveniently written as the system 



VPxx-\p 2 x + 2Vp 2 + \ = 0, [ZAl) 

where the second equation above (Ermakov-Pinney) defines V. A direct 
consequence of Ermakov-Pinney is 

B ■ p = 0, B = R-d x = d x + AVd x + 2V X - (2.12) 

For a summary of results and references on the Ermakov-Pinney equation see 



O . Having obtained the equation (|2.9| ) from (2.1) by using the reciprocal 



transformation ([O]), one may check that the standard WTC Painleve test 



[22 1 is satisfied. In fact the original equation ( |2.1| ) admits expansions with 
square root branching, but such generalized (weak |L8|| ) Painleve expansions 
were specifically excluded from the classification of Further details of 
the Painleve analysis will be given elsewhere |H|. 



3 The Lax pair 



Equation (|2.9|) , which is equivalent to the system ( 2.11 ), is integrable since 
it admits the Lax pair 

^xxx + Wifrx + (2V X - A)V = 0, 

IpT + ^~ 1 {p 2 1pXX -PPXtpX ~ {PPXX -P X + \)i>) = 0- 

The third order operator above is the Lax operator of the Kaup-Kupershmidt 
(KK) hierarchy [|ll|] , and so the system ( |2.11| ) may be considered as the first 
negative flow in this hierarchy. As far as we are aware this system has not 
been studied before, although it turns out to be an integrable generalization 
of the Tzitzeica equation. (The constant 2/3 in the T part of ( |3,ip can be 
gauged away, but allows for easy comparison with reference JO])- 

However, it is not the most general negative KK flow derived from such 
a Lax pair with an inverse power of A in the T part. Indeed, starting from 
a more general time evolution of the form 

ih + A" 1 (w-tPxx - \w x ^x + \{W X x + lQVW)^j = 0, (3.2) 
which is a reduction of the time part of a 2+1-dimensional non-isospectral 



KK hierarchy introduced in [11], the compatibility conditions of the third 



order Lax equation with (pj) yi e Ws the two conditions 



V T = ^W x , B ■ W = 0, (3.3) 
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where B = d x + 6 (d x Vd\ + d\ Vd x ) + 4 {d\ V + V d\ ) + 32 (d x V 2 + V 2 d x ) ■ 
The recursion operator of the KK hierarchy is 



Rkk 



Bd^Bdx 1 , 



(3.4) 



with B given by ( 2.12j ); B is a Hamiltonian operator for the KdV, KK and 
Sawada-Kotera hierarchies, (see [11] and references). 

Thus we see that the more general inverse KK flow ( p.3j ) is a system 
of higher order than ( 2.1 1] ) (while the most general first negative flow is 
Rkk ■ Vt = 0, of even higher order). However, the striking thing is that 
system ( 2.11| ) (or equivalently is a consistent reduction of (3J3). This 

is a consequence of the following remarkable operator identity: 



B ■ p 2 = 2(pd 2 x + bp x d x + lOpxx + 16V» ■ B ■ p. 
Thus for an arbitrary function h(T) we may substitute 



W 



p 



V 



Pxx Vji_ 
2p Ap 2 



KT) 
4p2 



(3.5) 



(3.6) 



into the T part (|3.2j ) of the Lax pair, and then upon scaling h(T) = 1 
(always possible for h ^ 0) the compatibility condition of the Lax pair 
becomes (|2.11 ). This is consistent with the second condition of ( |3.3[) due to 
the identity ( |3.5| ) together with B-p = ( |2.12| ). In the special case h(T) = 0, 
which is not relevant to the reciprocal transformation of the equation ( |2.1| ), 
we can integrate the resulting system to obtain 



p(\ogp) XT + p 3 + h(T) = 



(3.7) 



with h arbitrary. For h ^ we rescale and redefine T so that h = 1, in 



which case (3.7) it is known as the Tzitzeica equation [21], or with dependent 
variable <fi = logp the Bullough-Dodd equation [p| 



iiT + e 2 ^ + e" 



0. 



(3.8) 



If h = we have Liouville's equation. 

Finally we present the Lax pair for the original equation ( |2.lD (or equiv- 



alently (1.6)). Setting 



* = p 



(3.9) 



recalling that m 
that satisfies 



-p 3 , and transforming all derivatives in ([O), we find 



- ^xxx ~ '"A* 



{U X + ^) * 



0, 
0. 



(3.10) 



(As before the final constant 2/(3A) can be removed by a gauge transforma- 
tion.) 
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It is worth considering the effect of a Galilean boost and rescaling on the 
equation (|1.6| ), which restores some of the free parameters in (|1.4| ), yielding 

mt + cqu x + ju xxx + m x u + 3mu x = 0, m = u — a 2 u xx . (3-11) 

In that case, with a rescaled spectral parameter p = a 2 \, the Lax pair is 

(l-a 2 ^)*, = M (m + f + 

*t + J^x + (u - ^\ V x - U X ^ = 0. 

In the dispersionless case, if cq = = 7 then up to scaling the equation 
( 3.11| ) is equivalent to the Riemann shock wave equation 

ut + uu x = (3.13) 

in both the limits a 2 — > and a 2 — > 00, i.e. in the limit of both large and 
small wavenumbers.Q 

The generalized Tzitzeica equation has a variety of interesting so- 
lutions, including solitons on a constant background and a reduction to a 
Painleve III transcendent. We will treat these solutions in more detail in ||]. 
Here we mention that the reciprocal transformation takes the one-soliton so- 
lution on constant background to a soliton solution of the equation ( |3.11| ) 
with dispersion and vanishing at infinity, while the peakon solution ( |1.2| ) 
can be obtained from the soliton in the limit of zero dispersion (although 
the transformation Q2.7| ) breaks down for the peakon, when m = —p 5 is a 
Dirac delta function). 



4 Conservation laws 

In deriving infinitely many conservation laws for the equation ( |1.6| ) , we first 
introduce the quantity 

p=(logV)«- (4-1) 
Then from the spatial part of the Lax pair ( 3,1C| ) we have 

Cp := (1 - d 2 x )p = 3p Px + p 3 + Am, (4.2) 

while (removing the 2/(3A) term) the time part yields a conservation law 
for p, namely 



Pt=Jx, J=u x 



up-\- 1 {p x + p 2 ). (4.3) 

The density p, given by ( |4.1| ), may be written as a formal series in powers 
of the spectral parameter A, with coefficients determined recursively from 



1 Actually, the limit a 2 — * 00 of (3.11) in the dispersionless case is d 2 (u t + uu x ) = 0. 
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(4.2). Substituting a corresponding expansion for the current j into ( |4.3| ) 
and comparing powers of A yields an infinite sequence of conservation laws. 
However, it turns out that two different expansions are possible, leading to 
two infinite sequences of conserved quantities for JOf ). 

With m = — p 3 as before, the first expansion takes the form 



n=0 

with the corresponding expansion of the current being 



J 



n=0 



Thus the leading order (£ 1 term) in (p~3|) is just the conservation law (|2.6|). 



The first two densities and currents found recursively from (4.2, 4.3 ) are 

p {0) = -PxP' 1 , p {1) = IpxxP' 2 - pIp' 3 + 3^; UR . 

j(0) =Ux _ up (0) } = _ p 2 _ up (l) _ ' 

For brevity we omit the recursion relations for p^ n \j^ n \ but note that the 
densities take the form p^ n > = P^ /p 2 ™+ 1 with p( n ) a polyonmial of degree 
n + 1 in p,p x , . . . ,P( n +i)x', the even terms p( 2n ) are exact x derivatives. The 
densities p and p^ correspond to the conserved quantities and in 
(5.1) below. 

The second expansion is in positive powers of A, viz 



p = ^ r (n) A n+l ) J= ^ s Wr+l (4.7) 

n=0 n=0 

Again details of the recursive formulae for r^ n \g^ are omitted, but we note 
that the odd densities r^ 2n+1 ^ are all exact derivatives, and apart from 

r (0) = u> r (2) = u 3 ( 48 j 

(which yield the conserved quantities Hq and iT_i below) all other non- 
trivial densities in this sequence are non-local since they arise by repeated 
application of the inverse of the Helmholtz operator C = (1 — d 2 ) appearing 
on the left hand side of (fO|). 

There is another conservation law which does not appear in the above 
sequences, namely 

[v 2 xx + hv 2 x + 4v 2 ^ = (au 2 v - ^u 3 - AvC-\u 2 ) - v x £- l (u 2 ) x ^ , (4.9) 

where v is defined by 

v-iA-dD^u. (4.10) 

This conservation law yields the conserved quantity H\ in (|5.1| ) below, since 
the density on the left hand side of (|4.9| ) differs from the compact expression 
mv by a total derivative. 
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5 Hamiltonian structures 



As derived above, the equation fll.q ) has an infinite sequence of conservation 
laws. Here we list some of the simplest associated conserved quantities: 

= — | J* u 3 dx, Ho = — | J mdx, 

Hi = Ifmvdx, H 5 = J m 1 / 3 dx, (5.1) 

H 7 = -\ f(mlm~ 7/3 + 9m~ 1/3 ) dx. 

In the above we take u = (1 — d^)~ 1 m as before, and from ( f4.10|) we have 
v = (4 — — e?;ji) m. The labelling is such that H^ generates a flow 

of weight k with Hamiltonian operator B$ as in ( |5.2|) below. 

If we introduce the skew-symmetric differential operators Bq,B\ accord- 
ing to 

B = d x (l-d 2 x ){A-d 2 x ), B 1 = m^dxm^idx-d^m^dxm 2 / 3 , (5.2) 



then we can immediately write the equation ( |1.6| ) in terms of the gradient 
of a conserved quantity in two different ways, namely: 

SH , . 

771^50—^ = ^1—^. 5.3 



The identity (5.3) is our proposed bi-Hamiltonian form of the equation. 



However, in order to assert that the pair Bq, B\ define a genuine bi-Hamiltonian 
structure it is necessary to check that both are Hamiltonian (Poisson) oper- 
ators, and that they are compatible [17, §. The first operator Bq is constant 



coefficient and so the Jacobi identity is trivial and it is clearly Poisson, but 
for the non-local operator B\ we have not succeeded in verifying the Jacobi 
identity. 

We have noted other properties of the operator pair ( |5.2[ ) which strongly 
suggest that they do define the correct bi-Hamiltonian structure for (|1.6| ). 
In the list ( |5,ip we see that Hq is a Casimir for Bq and H§ is a Casimir for 
B\. The translational flow in the hierarchy is 

5Hi 5H 7 , . 

m x = B — ^ = B 1 —L, 5.4 
dm dm 



where the second equality in ( |5.4| ) requires suitable definition of the integral 
operator d~ x in B\. Also, the fifth-order vector field Bq \ s j us ^ ^he recip- 
rocal transformation of the fifth-order Kaup-Kupershmidt equation. This is 
strong evidence that the operator Bq B^ 1 should be the reciprocal transform 



of the Kaup-Kupershmidt recursion operator (3.4). 



9 



6 Peakon dynamics 



Here we consider the dynamics of multi-peakon solutions of the equation 



(1.6), which take the form 

N N 

u(x,t) = Y,Pj( t ) e ~ lx ~ qm , m = 2j2p j (t)8(x-q j (t)), (6.1) 

3=1 3=1 

with N, the number of peakons, being arbitrary. 

Substituting the expression fl6,l|) into the equation ( |1.6[) yields the cou- 
pled ODEs for qj(t),Pj(t), namely 

Pj = 2 J2k=iPjPksgn( qj - q k )e-\v- qk \ . , 



k=iPke 



which is just the system ( |1.8| ) in the particular case 6 = 3. Although the 
variables qj,Pj are not canonical, and we have not found a (non-canonical) 
Hamiltonian structure except in the special case N = 2, the dynamical 
system (6.2) is integrable in the sense that it arises from a Lax equation 



L=[M,L] (6.3) 

where L and M take the form 

L = I(A-A T -1-2C)P, M = D-(A-A T )P. (6.4) 

In order to specify the matrices in the Lax pair, which are built out of 
the strictly lower triangular matrices A, C and the diagonal matrices P, D, 
we first note that an immediate consequence of the system fl6.2|) is that 



^sgn(qj — q k ) = 0, so that the peakons preserve their relative ordering 
(as for Camassa-Holm peakons). Without loss of generality we take the 
ordering q\ < t/2 < • • • < <7iv and then for the matrices appearing in Q6.4| ), 
the non-zero components of A and C may be given (without modulus signs) 
as 

A jk = e«*-V, C jk = l, j>k, (6.5) 
while the diagonal elements of P, D are the components of the vectors 

p = ( Pl ,p 2 ,..., PN ) T , d = (A-A r )p (6.6) 

respectively. 

In the case of 2-peakon scattering, i.e. the equations ( |6.2[) for N=2, we 
consider two peakons that are initially well separated with q\ < qi and have 
asymptotic speeds c\ and C2 as t — > — oo with c\ > C2 and c% > so that 
they eventually collide. The situation is shown in Figure ffl. 
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u(x,t) u(x,t) 




Figure 1: Two rear-end collisionss of b = 3 peakons. The initial positions 
are x = 33 and x = 66. The faster peakon moves at three times the speed 
of the slower one. For this ratio of speeds, both collisions result in a phase 
shift to the right for the faster space-time trajectory, but no phase shift for 
the slower one. 



In this case, we have two conserved quantities 

P ■= P1+P2 = C!+c 2 , H := pl+pl+2p x p 2 e qi - q2 {2-e qi - q2 ) = c\+c 2 2 . (6.7) 

The equations of motion are then most easily solved in terms of the sum 
and difference variables P, Q = q\ + q 2 , p = Pi — P2 arid q = q\ — q 2 , and 
the quadratures involved are almost the same as in the Camassa-Holm case 
(see For the positions of the peaks we find the explicit formulae 

qi(t) = cit + ilog(4r(ci-c 2 ) 2 )-log(r e ( Cl - C2 )' + 4 C2j p), (6.8) 
q 2 (t) = C2 t-ilog(4r(ci-c 2 ) 2 )+log(re( Cl - C2 )'+4 ClJ p), (6.9) 

where T is an integration constant. 

In the limit t — > oo, these formulas show that the solitons exchange their 
asymptotic speeds, or equivalently their momenta and amplitudes, as 

lim qi(t) = c 2 t and lim q 2 {t) = c\t . (6.10) 

t—*oo t^oo 

Thus, the main effect in the peakon scattering is an exchange of momentum, 
or amplitude, between the two peakons, resulting only in a phase shift at 
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asymptotic times. The phase shift for the faster soliton (the one with speed 
c\ in the limit t — > -co) is defined and evaluated using (|6.8[) and (|6.9|) as 



A?/ = <?2(+oo) - gi(-oo) = log 



cijci +c 2 ) 
(ci - c 2 ) 2 



(6.11) 



Likewise, the phase shift for the slower soliton (the one with speed c 2 in the 
limit t — > — oo) is defined and evaluated using (rcTq) and (^313) as 



A<7s = gi (+00) - 92 (-00) = log 



(ci - c 2 ) 2 
c 2 (ci + c 2 ) 



(6.12) 



So for ci/c 2 > 3 both peakons experience a forward shift, and for 1 < ci/c 2 < 
3 the faster peakon shifts forward while the slower one shifts backward; the 
case ci/c 2 = 3 is the turning point at which the slower peakon has no phase 
shift. Remarkably, phase shift scattering rules corresponding precisely to 
these hold for Camassa-Holm peakons, except with 3 — > 2 ||]. 

An interesting special case is the peakon-antipeakon collision when c\ 
and c 2 have opposite signs; only in this case can the solitons overlap when 
the variable p = Pi — P2 diverges. For the perfectly antisymmetric collision 
ci = — c 2 = c the resulting solution of the partial differential equation ( 2.1 
is 



u(x, t) - 
m(x, t) 



2c 



-2c\t\ 



e -\x+c\t\ \ _ 

5{x + c\t\) 



\x— c\t\ 



S(x - c\t\) 



(6.13) 



In numerical simulations, we also investigated the emergence of peakons 
from a Gaussian initial condition. The integrable behavior is evidenced 
in Figure |2| as the peakons collide elastically as they recross the periodic 
domain. 



7 Conclusions 

In a forthcoming article || we will consider an integro-differential generaliza- 
tion of (|Q|), extending the results of pi to include an extra parameter 6, such 
that there are multi-pulson solutions described by the finite-dimensional sys- 
tem ( |1.8[) with the generating function 

1 N 

Gn = 2 ^2 Pip i 9 ( Qi ~ q ^ ■ 

for an arbitrary even function g[x) (so that g(x) = e _ ' x ' in the peakons case). 
It turns out that for arbitrary b values the two-pulson (N = 2) dynamics is 
Hamiltonian and integrable. 
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u(x,t) u(x,t) 




10 20 30 40 50 60 70 80 90 100 



X 

Figure 2: b = 3 peakons emerging from a Gaussian of unit area and a = 5 
centered about x = 33 on a periodic domain of length L = 100. The fastest 
peakon recrosses the domain five times and has many elastic interactions 
with the slower ones. 
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